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Abstract 

Some indication of conditions that are necessary for the formation of black 
holes from the collision of bubbles during a supercooled phase transition in 
the the early universe are explored. Two colliding bubbles can never form a 
black hole. Three colliding bubbles can refocus the energy in their walls to 
the extent that it becomes infinite. 
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I. INTRODUCTION 



A small number of black holes being produced in the early stages of universe would have 
important cosmological consequences, both from their contribution to the density of the 
universe and to the high energy cosmic ray background For this to have occured there 

would have to have been density perturbations with large amplides because small amplitude 
perturbations do not grow in size in the primordial plasma. The isotropy of the universe 
suggests to us that the universe was fairly homogeneous on the length scales corresponding 
to galaxies or larger objects, but comparatively little is known about the homogeneity of 
the early universe on small scales, for example scales much less than the horizon scale at 
nucleosynthesis which encompassed a mass of around 3 x 10 36 Kg. 

One situation in which a substantial amount of inhomogeneity may have arisen could be 
at a first order phase transitions where a gauge symmetry was broken. The possibility of 
black hole formation was raised very early in the theory of these phase transitions [f|,0. 

First order phasel) transitions are characterised by tunnelling from the old to the new 
phase. For spontaneous symmetry breaking the phases are distinguished by the values of a 
Higgs field. The tunnelling is modeled by the instantaneous nucleation and expansion of a 
bubble of the broken-symmetry phase |6|j7|]. In the supercooled case, where the tunnelling 
is smaller than the expansion rate of the universe, it is possible for the thermal energy of 
the old phase to be surpassed by the energy of the vacuum, then the bubble walls accelerate 
to speeds close to the speed of light and build up considerable amounts of energy. 

The collision of these bubbles and the possibility of forming black holes was considered 
in reference f|. We found that when two bubbles collide the bubble walls could continue 
through one another unchanged or the energy could convert into the phase of the Higgs 
field. In the latter case, most of the energy emerged from the collision region as phase waves 
moving at the speed of light. We also argued on energy grounds that eight colliding bubbles 
could form a black hole. 

It is possible to model the collision of two thin-walled bubbles by distributional gravi- 
tational sources ||. The bubbles start out with flat Minkowski inside and de Sitter space 
outside. Each bubble is invariant under Lorentz boosts and rotations |9[], therefore the two- 
bubbles with a prefered axis have the symmetry 0(2,1) ||. This has as many killing fields 
as spherical symmetry and Wu was able to demonstrate that there is a modified form of 
Birkoff 's theorem that can be used to find the metric. 

The rapid acceleration of the bubble walls up to luminal speeds indicates that the collision 
can be simplified by replacing the bubble walls with null surfaces. When the energy from 
the walls is converted into phase waves, as often happens, the energy also leaves the collision 
on null surfaces. The collision of two bubbles becomes the collision of two null shock waves, 
a problem whose solution is known ||10Ull|| . 

By making these reasonable simplifications it becomes possible to include a third or a 
fourth bubble in the collision. Under suitable conditions we shall see that the energy of the 
bubble walls is sufficient to refocus itself into a singularity and therefore form a black hole 
or a naked singularity. 
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II. 0(2,1) INVARIANT METRICS 



The group 0(2,1) is the symmetry group of the pseudosphere H 2 - We can obtain a 
spacetime metric which is invariant under 0(2,1) by replacing the sphere by the pseudosphere 
in the usual 0(3) invariant metric, 

g = -S(s, x)ds 2 + X(s, x)dx 2 + s 2 (d6 2 + sinh 2 Odcj) 2 ) . (2.1) 

Substituting this form of the metric into the vacuum Einstein equations (with a cosmological 
constant 3/a 2 ), we can deduce easily that X — 1/S — f, where 

2M s 2 j , 

/ = 1 + -• 2.2 

s a 2 

The parameter M is a constant but it does not have the same physical associations of a 
mass. More details about these metrics can be found in reference ||12j| . 

We can draw Penrose diagrams of the surfaces that are orthogonal to the pseudspheres. 
The case with vanishing cosmological constant, which will be called pseudo-Schwartzchild, 
is shown in figure |l|. This is the usual Penrose diagram of Schwartzchild rotated through 
90°. Surfaces of constant coordinate s are spacelike in s > 2M and the metric aproaches 
flat space (Milne universe form) in the limit that s —>■ oo. 

The metric with M = is de Sitter space. It will be useful to have the embedding of the 
metric into a five dimensional Minkovski spacetime with cordinates X = (X, Y, Z, W, V), 
where de Sitter space is defined by the hyperboloid 

X ■ X = X 2 + Y 2 + Z 2 + W 2 - V 2 = a 2 . (2.3) 

The two de Sitter coordinate systems are related by 

X = s sinh 6 cos 0, (2.4) 

Y = s sinh 6 sin <ft, (2-5) 
Z = (a 2 + s 2 ) 1/2 sm(x/a), (2.6) 

W = (a 2 + s 2 ) 1/2 cos(x/a), (2.7) 

V = scosh9. (2.8) 

During a supercooled phase transition the universe becomes locally similar to de Sitter 
space. Bubbles of the true vacuum phase (Minkowski space) appear instantaneously. Their 
walls move along trajectories with 0(3,1) symmetry, 

(X - X ) • (X - X ) = r 2 . (2.9) 

The radius r$ is the initial bubble radius and is fixed by physical parameters. For phase 
transitions at temperatures far below the Planck temperature of 10 19 Ge^, this radius is far 
smaller than the cosmological horizon size a. 

If ro = we can replace the bubble surface by a null surface. The generators of this 
surface are light rays 

X-X = AP, (2.10) 



3 



where A is an affine distance along the ray and the constant momenta P satisfy P • P = 0. 
This is the approximation that will be used henceforth. The light rays can be also be 
parameterised in an 0(2, 1) invariant form by the substitution X(s,x, 9, (ft). Differentiating 
with respect to the affine parameter gives 

s = s- 1 (ZP Z + WP W ) , (2.11) 
x = (sf)' 1 (WP Z - ZP W ) , (2.12) 
<P = s~ 2 csh 2 # (XP Y - YP X ) , (2. 13) 

where the components of P are denoted by (P x , P Y , P z , P w , P v ). The values of s etc. are 
the components of the tangent vectors in the 0(2, 1) coordinate frame, with the latter two 
quantities in brackets being the conserved linear and angular momenta respectively. 



III. TWO BUBBLE COLLISIONS 

In this section we will take a look at the problem of two intersecting null shells, being 
relevant to the collision of two bubbles when the energy from the bubble walls leaves the 
collision region at the speed of light. Numerical studies have shown that this is a good 
approximation when there is a large difference of phase of the Higgs field in each of the 
bubbles. 

Consider a null hypersurface with null geodesic generator 1 and affine parameter u. It 
is convenient to introduce a 'pseudonormal' null vector n = d/dv, normalised by l.n = — 1 
and parallel propagated along 1. The remaining tangent vectors ei, % — 1,2, define a two 
dimensional surface with induced metric 

Kb = 9ab + kn b + n a l b (3.1) 

and it is convenient to have them parallel propagated along 

The embeddings of the null surface are described by extrinsic curvature tensors k ab and 
m a b, where 

k ab = h a c h b d n c . d (3.2) 
m ab = h a c h b % d . (3.3) 

From the condition that e ; • 1 = e ; • n = 0, it is also possible to write 

ka^ef = -n.(Viej), (3.4) 
= -l.(Viej) (3.5) 

Their contractions will be denoted by k and m. It is easy to show that m is also the 
expansion of the null geodesic congruence, that is m — V • 1. 

A distributional source of stress-energy on the shell would take the form 

Tab = <JlJbS(v). (3.6) 
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Junction conditions across shell can be obtained from integrating the Einstein field equations. 
(See ref [13| for example, but note that 1 and n have the opposite meaning.) These force the 
metric to be continuous and allow the choice of 1, n and ei to be continuous. Furthermore, 
they relate the extrinsic curvatures of the surfaces to the stress energy. If the surface is 
moving to the right from a region C into another region B then, 

87rGcr = ks — kc, and m# = mc- (3.7) 

where G is Newton's constant. 

Two null hypersurfaces divide spacetime into four regions, labeled as in the diagram |2|. 
The null geodesic generators of the surfaces % = 1 ... 4 are labeled by li and the extrinsic 
curvatures will be indexed accordingly. In the vicinity of the intersection pseudosphere, the 
junction conditions imply that 

SnG^ = k l c - k l D , 8nGa 2 = k 2 c - k% (3.8) 
8nGa 3 = k%- k 3 A , 8nGa 4 = k 4 D - k\ (3.9) 

It is also possible to obtain normal vectors from the Z's, for example 

ni = -(li • la) -1 l a and n 2 = -(U • l 2 ) %■ (3.10) 

These relations are specific to the regions in which the vectors are defined, shown in fig. |2|, 
but the normal vectors rii are themselves continuous across the i'th surface. The continuity 
equations for the m's then imply the following equations for the fc's, 

(h ■ l a )A& = (l a • la)*!, (U " h)kl = (h ■ h)k% (3.11) 
(la • l 2 )k% = (la • U)k% (li • h)k 2 c = (li • h)kj>. (3.12) 

Finally, we allow the affine parameterisations on either side of the collision region to differ 
by a constant scaling, 

1 3 = 7I1, h = (3h- (3.13) 

The continuity equations can be written in a much more useful form. Introduce 

Fc = (h-h)khk 2 c , (3.14) 

using extrinsic curvatures of surfaces adjacent to region C, and similarly for each of the 
regions in turn. Then, 

F C F A = (h ■ 1 2 )(1 3 • U)k l c k 2 c k\k A Al (3.15) 
Using the continuity relations gives, 

F C F A = (h ■ 1 4 )(1 2 • iMkjjklkl = F B F D . (3.16) 



This will be refered to as 'the pseudo-DTR relation' fllQlJTTf . 

The metrics that are of interest to us have the 0(2,1) symmetry. Continuity of the 
pseudosphere sections across the null hypersurfaces implies that the 9 and 4> coordinates are 
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continuous, and also the s coordinate because it sets the area element, but the x coordinate 
will be different on either side of the surfaces. 
The null generators of the surfaces are then 

h = e s -r 1 e*, I 2 = e B + /- 1 e x . (3.17) 

Using the connection components of the metric gives, 

-k 1 c = k 2 c = f c /s. (3.18) 

Therefore the pseudo-DTR relation reads 

fcfA = f B fo. (3-19) 

This is the junction condition that determines the metric in any region if the metric is known 
in the other three. It is also possible to get the surface energy density from 



a 3 = ^— f f (3.20) 

For the collision of two bubbles, the initial region A is de Sitter space, 

f a = 1 + s 2 /a 2 . (3.21) 

Regions B and D are inside the bubbles, where the metric is flat. The energy on the bubble 
surface follows from equation |3.20| . Noting that the affine length along the surface is A = s, 
this can be written 

a = (SitGa 2 )-^. (3.22) 

The increase in surface energy represents the latent energy of the false vacuum that the 
bubble absorbes as it expands. 

The intersection region C is in the true vacuum where there is no cosmological constant 
and the metric is pseudo-Schwartzchild, 

f c = 1 - 2M/s. (3.23) 

The collision takes place at a particular value of s, say s c . From the pseudo-DTR relation 
we have 

(l + s 2 Ja 2 ) (1 -2M/s c ) = 1. (3.24) 

This gives 

2M = -5T-5- ( 3 - 25 ) 

In particular, we see that M is always positive and we always get s c > 2M. This is the 
condition for there to be no singularity anywhere in the spacetime. Therefore, in this 
simplified form of the problem, the collision of two bubbles never produces a black hole. 
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IV. MANY BUBBLE COLLISION 



Although the collision of two bubbles fails to produce a black hole, the collision can 
come quite close to the limit s = 2M where the spacetime becomes singular. Only a small 
additional amount of energy from another bubble should be required to tip the balance 
in favour of a black hole. Unfortunately, introducing a third bubble breaks the symmetry 
and prevents us writing down the metric explicitly. What we can do instead is aim to find 
conditions under which singularities are formed. 

As a first step we shall consider the trajectories of light rays in the spacetime formed 
from two colliding bubbles, especially their refraction from the bubble walls (fig. We 
shall write the momentum vector p of a general ray in terms of the generators of one of the 
bubble walls, 

p = p u \ + p v n + p 1 e 1 +p 2 e 2 . (4.1) 

If the light ray is itself part of a bundle from the same initial point, then the expansion 9 of 
the bundle satisfies Raychaudhuri's equation, 

p(0) = -R ab p a p b - 2a 2 - \d\ (4.2) 

where a is the shear and p has been identified with its directional derivative following a 
common practice. (If we define the analogue of m ao for the congruence generated by p, 
then 9 and a are the trace and trace free parts of m a &.) The Einstein equations relate the 
curvature to the surface energy in the bubble wall, 

T ab = (rUbS{v), (4.3) 

but we also have 

89 89 

Integrating gives the focusing rule, 

9 B -9 A = 8nGal-p, (4.5) 

for a ray passing from region A into region B. 

We also see from Raychudhuri's equation that the expansion is finite, suggesting that 
the momentum vector should be continuous. Nevertheless, there is some refraction of the 
ray in the 0(2,1) coordinate system, due to the discontinuity of the coordinates themselves. 
Consider again the ray that pases from region A to region B, but now using the explicit 



forms of the bubble generators |3.17|, 



P = |(s - fAx)h + |(s + f A x)U + 9e e + 0e . (4.6) 
In region B we use the vector 1 2 , related to I4 by equation |3.13j , 



U = (f B /fA)h. (4.7) 
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Restoring the s and x coordinates gives 



s \ i (l+P (3- 
fx] B H/?-l l + /3J\fx 

where (3 = /b//a- Similarly, for the passage from B to C, 

s \ i ( 1 + 7 1 — 7 \ / s 



(4.8) 



/x7 c Hi-7 i + -yJ\fiJ B ' (4 ' 9) 

where 7 = / c //b- 

A particular case of interest for us later on is where the ray passes close to the collision 
region, and the values of s in both of the refraction relations are the same. Then the 
combination gives 



fx) c \0 (3j\fx 



(4.10) 



§ c = e A - ^Gaf A l (s 2 - fx 2 ) A . (4.11) 



and 



Now we turn to the collision of three bubbles in a symmetrical arrangement with centres 
equally spaced a proper distace d appart, as shown in figure |j. The cartesian coordinates 
can be chosen so that the bubbles all begin at the same value the time coordinate V and in 
the same plane. We therefore take the centres to be at, 

X 1 = (0,0,^,^1,0), (4.12) 

X 2 = (0,0,-^,^,0), (4.13) 

X 3 = (X 3 , 0,0, W 3 ,0), (4.14) 

where 

Z x = asin(d/2a) (4.15) 

and 

XI = (4,6, 

We also have the condition that Xj • X, = a 2 to fix W± and W3. 

Bubbles 1 and 2 nucleate at s = in the 0(2,1) coordinate system. The bubble surfaces 
are given by the ray equation |2.10| and collide on a surface of constant coordinate s, 

^ (4-17) 



a 2 - Z 2 



Some elementary geometry shows that the generators of bubble 3 intersect bubble 2 after 
an affine length 
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Substituting for X 3 gives 



2Z 2 

Xp= WiPW + Z^z- (419) 



This also gives a value of the coordinate s for the intersection, 

p ~ Wt W 1 P W + Z 1 P Z ' 1 } 

We shall restrict attention to the symmetrical rays that pass through the collision points. 
For these rays, x = P z = 0, and from equation |2.11| , 



pW 

s P = ^- (« 2 + Zl) . (4.21) 



The surface energy of bubble 3 at the intersection points is given by equation |3.22 



(4.22) 



0a = - = (4-23) 



p ^GaW x P w ' 
The expansion of the bubble wall is given by 

n W X P W 

X p z\ 

Now it is possible to use the focussing rule [4.11| , 

6c = 6a ~ SnGapfX 1 ^ (4.24) 

This gives, 

6 c = ^^{a e -2Zl(a 2 + Z!) 2 }. (4.25) 
Therefore the expansion of the outgoing bubble wall is negative after the collision when 



Zi > za, where z « 0.54512 is a root of z 3 + z - Jl/2 = 0. If Z\ > ay/3/2, then the de 
Sitter geometry does not allow the three bubbles to have a common collision point and there 
is no refo cussing. 

Raychaudhuri's equation tells us that once the expansion becomes negative it remains 
negative and diverges in a finite affine distance along the ray. When the expansion diverges 
this forces the surface energy to become infinite. At this point there is a spacetime singular- 
ity, in the sense that the curvature is more singular than the distributional curvature with 
which we began. 

The preceeding arguments have established that refocussing occurs along a line extending 
from the collision point in the direction of the symmetrical ray. This line was chosen only 
for convenience and the refraction rules obtained earlier can be used to examine refocussing 
anywhere on the bubble wall. 
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V. CONCLUSION 



We have argued that the bubbles nucleating at a supercooled phase transition can be 
replaced by expanding null shells of energy. With this simplification, the collision of two 
bubbles never produces a singularity, in agreement with previous results. 

When three bubbles collide the surface energy in parts of the walls can be refocused to 
the extent that it becomes infinite. What happens next depends on general issues of general 
relativity such as the cosmic censorship conjecture. If the singularities are to be invisible 
from future infinity of the emerging universe then they would have to be hidden behind an 
event horizon, and form a black hole. If cosmic censorship fails, then the bubble collision 
results in a new object, but one whose time-evolution presents problems in general relativity. 

The condition for refocusing the energy along a symmetry axis for a symmetrical collison 
was a proper separation d > 1.153a. These three bubbles have a triple collision point 
provided that d < 2.094a. Although the non-symmetrical case is more complicated it is 
theoretically possible to follow it through using the equations that have been given. In 
particular, small deviations from the symmetrical case should leave the results unchanged. 

If there are four colliding bubbles, then there is a question of how close the boundary 
of the intersection region E in figure f| is to being future trapped. The arguments given 
for three bubbles can be generalised to show that the symmetrical null generators of this 
surface are converging for large bubble separations. Of course, the surface will never be 
converging at the (one dimensional) corners because of lack of smoothness. This raises two 
interesting questions. If the surface where trapped, then would it be possible to generalise 
the singularity theorems to conclude that the singularity is generic, given that the strong 
energy condition only holds to the future of the first collision? What is the effect of replacing 
a smooth trapped surface with a piecewise smooth surface? 
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FIGURES 



FIG. 1. Penrose diagram of the pseudo-Schwartzchild spacetime with M > 0. The s coordinate 
increases from bottom to top. 

FIG. 2. Two intersecting null hypersurfaces. 

FIG. 3. Refraction of a light ray with momentum p by two bubbles. 

FIG. 4. The collision of three bubbles and a ray from X3 that passes through the bubbles 
symmetrically. 
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